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ABSTRACT 

We consider the global cosmological evolution and the evolution of the density con- 
trast in the Bose-Einstein condensate dark matter model, in the framework of a Post- 
Newtonian cosmological approach. In the Bose-Einstein model, dark matter can be 
described as a non-relativistic, Newtonian gravitational condensate, whose density and 
pressure are related by a barotropic equation of state. For a condensate with quartic 
non-linearity, the equation of state is polytropic with index n = l.Thc basic equation 
describing the evolution of the perturbations of the Bose-Einstein condensate is ob- 
tained, and its solution is studied by using both analytical and numerical methods. 
The global cosmological evolution as well as the evolution of the perturbations of the 
condensate dark matter shows significant differences with respect to the pressureless 
dark matter model, considered in the framework of standard cosmology. Therefore the 
presence of condensate dark matter could have modified drastically the cosmological 
evolution of the early universe, as well as the large scale structure formation process. 

Key words: cosmology: theory: dark matter: large-scale structure of Universe - 
instabilities-equation of state 



^ 1 1 INTRODUCTION 



Cosmological observations provide compelling evidence 
, that about 95% of the content of the Universe resides in 
■ two unknown forms of energy that we call dark matter and 
dark energy: the first residing in bo und objects as non- 
luminous matter JPersic et al.l 1 19961 : iBoriello and Saluccil 
l200ll : iBinnev fc Tremaind 120081 ). the latter in the form 
of a zero-point energy t hat pervades the wh ole Universe 
jPeebles fc Ratral 120031 ; IPadmanabhanl 120031 '). The dark 
matter is thought to be composed of cold neutral weakly 
interacting massive particles, beyond those existing in 
the Standard Model of Particle Physics, and not yet 
detected in accelerators or in dedicated direct and indirect 
searches. There are many possible candidates for dark 
matter, the most popular ones being the axions and 
the weakly interacting massive particles (WIMP) (for 
a r eview of the particle physi cs aspects of dark matter 
see lOverduin fc Wesson! l|2004 )). Their interaction cross 
section with normal baryonic matter, while extremely 
small, are expected to be non-zero and we may expect 
to detect them directly. It has also been suggested that 
the dark matter in the Universe might be composed of 
superheavy particles, with mass > 10 10 GeV (|Chung et all 
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1 19981 . 120051 : iKolb et al.l [20071 : IChuzhov fc Kolbll2009l ). But 

observational results show the dark matter can be com- 
posed of superheavy particles only if these interact weakly 
with normal matter or if th eir mass is above 10 15 GeV 
jAlbuauerque fc Baudij l2003h . Scalar fields or other long 
range coherent fields coupled to gravity have also inten- 
sively been used to mode l galac ti c dark matter jLee fc Kobl 
199d: iNucamendi et ail |2000|; iMatos fc Guzmanl 12001 



Mielke fc Schunkl 120021: lArbev et all 120031: 



Fuchs 

~n — 



Mielke 



2004; Hernand ez et alj|2004 iGianniod 2005: Khlopo v et al 
20051: bernal fc Guzmanl |200d : lArbevH200rj , I2OO8I : lBrisc~ 



20111 ). Alternative theoretical models to explain th e galactic 



rotat i on curves have also been elaborated re c ently llMilgrom 
1 19831 : iMannheiml Il993l: [Bekensteinl 120041: iMak fc Harko 
12004 : iBrownstein fe Moffat] 120061: iHarko fc Cheng 



. iBertolami et all 120071; iBoehmer fc Harkol l20(m 
boehmer et al.ll2008al lbl: ICapozziello et al.ll2009l ). 

In order to explain the recent observational data, 
the ACDM (A Cold Dark Matter) model w as developed 
i Peebles fc Ratral 120031 : IPadmanabhanl [2003). The ACDM 
model successfully describes the accelerated expansion of the 
Universe, the observed temperature fluctuations in the cos- 
mic microwave background radiation, the large scale matter 
distribution, and the main aspects of the formation and the 
evolution of virialized cosmological objects. 

Despite these important achievements, at galactic scales 
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~ 10 kpc, the ACDM model meets with severe difficulties in 
explaining the observed distribution of the invisible matter 
around the luminous one. In fact, iV-body simulations, per- 
formed in this scenario, predict that bound halos surround- 
ing galaxies must have very characteristic density profiles 
that feature a well pr onounced central cu sp, pnfw{t) — 
p a /{r/r s )(l + r/r s ) 2 l|Navarro et alj ll997T ), where r a is a 
scale radius and p s is a characteristic density. On the ob- 
servational side, high-resolution rotation curves show, in- 
stead, that the actual distribution of dark matter is much 
shallower than the above, and it presents a nearly constan t 
density core: ps(r) = p rl/(r + r )(r 2 + rl) (iBurkertll 19951) . 
where ro is the core radius and po is the central density. 

At very low temperatures, all particles in a dilute Bose 
gas condense to the same quantum ground state, form- 
ing a Bose-Einstein Condensate (BEC), i.e., a sharp peak 
over a broader distribution in both coordinates and mo- 
mentum space. A coherent state develops when the parti- 
cle density is enough high, or the temperature is sufficiently 
low. The Bose-Einstein condensation process was first ob- 
served experimentally in 1995 in dilute alkali gases, such 
as vapors of rubidium and sodium, confined in a magnetic 
trap, and cooled to very low temperatures. A sharp peak in 
the velocity distribution was observed below a critical tem- 
perature, indicating that condensation has occurred, with 
the alkali atoms condensed in the same ground state and 
showing a narrow pe ak in the momentum space and in 
the coordinate space ([Anderson et al. I Il995l ; iBradlev et al 
1 19951 ; [5 avis et all ll995T ). Quantum degenerate gases have 
been created by a combination of laser and evaporative 
cooling techniques, opening several new lines of research, 
at the bord er of atomic, statis t ical and condensed mat- 
ter physics (iDalfovo et all Il999l; ICornell fc Wiemanl 1200$ 



Ketterld 120021; iPitaevskii fc Stringaril 120031; iDuine fc Stood 
2004IChen et alj|2005l ; iPethick fc Smitbll200gt ).~ 

The possibility that dark matter could be i n the form 
of a Bose-Einstein condensate was considered in 
and IJi fc Sir] <|l994 ). The condensate was described by the 
non-relativistic Gross-Pitaevskii equation, and its solution 
was obtain ed numerically. An a lterna tive approach was de- 
veloped in iBoehmer fc Harkol |2007t) . By introducing the 
Madelung representation of the wave function, the dynamics 
of the system can be formulated in terms of the continuity 
equation and of the hydrodynamic Euler equations. Hence 
dark matter can be described as a non-relativistic, Newto- 
nian Bose-Einstein gravitational condensate gas, whose den- 
sity and pressure are related by a barotropic equation of 
state. In the case of a condensate with quartic non-linearity, 
the equation of state is polytropic with index n = 1. To 
test the validity of the model the Newtonian tangential ve- 
locity equation of the model was fitted with a sample of 
rotation curves of low surface brightness and dwarf galax- 
ies, respectively. A very good agreement was found between 
the theoretical rotation curves and the observational data 
for the low surface brightness galaxies. Therefore dark mat- 
ter halos can be described as an assembly of light individ- 
ual bosons that acquire a repulsive interaction by occupying 
the same ground energy state. That prevents gravity from 
forming the central density cusps. The condensate parti- 
cle is light enough to naturally form condensates of very 
small masses that later may coalesce, forming the struc- 
tures of the Universe in a similar way than the hierarchical 



clustering of the bottom-up CDM picture. Then, at large 
scales, BEC perfectly mimic an ensemble of cold particles, 
while at small scales quantum mechanics drives the mass 
distribution. Different properties of the Bose-Einstein con- 
densate dark matter halos, like the effects of the rotation 
and of the vortices, as well as the cosmo logical effects of 
the condensation were also investigated (Ferrer fc Grifolj 
| 2004 iFukuvama et al] 120081; IFukuvama fc Morikawal \20(M 
Rindler-Daller fc Shapiro! 120091; ISikivie fc Yand 120091; I Lee I 



1 20091 ; iBrook fc Coles! 120091 ; iKain fc Lind |2010| ; iLee fc Lir 

boioh . 

It is the purpose of the present paper to study the 
global cosmological dynamics of gravitationally self-bound 
Bose-Einstein dark matter condensates, and the evolution 
of the small cosmological perturbations in the condensate. 
The equations of motion of the condensate dark matter are 
obtained in a Post-Newtonian approximation by using the 
conservation of the general relativistic energy-momentum 
tensor, and considering the small velocity limit. The cos- 
mological dynamics of the Bose-Einstein condensate is also 
studied. The exact solution of the Friedmann equations is 
obtained, and it is compared with the standard Einstein- 
de Sitter cosmological model. In order to study the evolu- 
tion of the small cosmological perturbations the equation 
describing the Newtonian perturbations with pressure is ob- 
tained in a general form, by also taking into account a term 
which was neg l ected in the previous studies of this problem 
(lMcCrealll95l; lHarrisonlll965l: iLima et all Il997l : iReisI 120031 ; 
lAbramo et~aT!fl2007l ; iPace et al.ll2010r ). The equation of the 
density contrast for the Bose-Einstein condensate is investi- 
gated by using both analytical and numerical methods. The 
presence of the condensate dark matter significantly mod- 
ifies the cosmological dynamics of the Universe, as well as 
the large scale structure formation. 

The present paper is organized as follows. The basic 
properties of the Bose-Einstein condensate dark matter ha- 
los are reviewed in Section [5] The Post-Newtonian hydro- 
dynamical equations of motion for a perfect fluid with pres- 
sure are derived in Section [3] The cosmological dynamics of 
the Bose-Einstein condensate dark matter is considered in 
Section [4] The equation describing the small cosmological 
perturbations of a fluid with pressure are derived in Section 
The evolution of the small cosmological perturbations in 
a Bose-Einstein condensate is considered in Section [6] We 
discuss and conclude our results in Section [7] 



2 DARK MATTER AS A BOSE-EINSTEIN 
CONDENSATE 

In a quantum system of N interacting condensed bosons 
most of the bosons lie in the same single-particle quantum 
state. For a system consisting of a large number of parti- 
cles, the calculation of the ground state of the system with 
the direct use of the Hamiltonian is impracticable, due to 
the high computational cost. However, the use of some ap- 
proximate methods can lead to a significant simplification 
of the formalism. One such approach is the mean field de- 
scription of the condensate, which is based on the idea of 
separating out the condensate contribution to the bosonic 
field operator. We also assume that in a medium composed 
of scalar particles with non-zero mass, when the medium 
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makes a transition to a Bose-Einstein condensed phase, the 
range of Van der Waals-type scalar mediated interactions 
among particles becomes infinite. 



2.1 The Gross-Pitaevskii equation 

The many-body Hamiltonian describing the interacting 
bosons confined by an external potential Vext is given, in 
the second quantization, by 



H 



dm + (r) 



2m 



V 2 + Vrot (r) + V ext (r) 



i J drdr $+ (r) <l + (f) V (r - r ) <£ if) $ (f) , 



(1) 



where $ (r) and $ + (f) are the boson field operators that 
annihilate and create a particle at the position r, respec- 
tively, and y (f — r 1 ) is the two-body interatomic potential 
[jPalfovo et al.ll 19991 : IChen et alJl2005h . V rot (r) is the poten- 
tial associated to the rotation of the condensate. 

The use of some approximate methods can lead to a sig- 
nificant simplification of the formalism. One such approach 
is the mean field description of the condensate, which is 
based on the idea of separating out the condensate contribu- 
tion to the bosonic field operator. For a uniform gas in a vol- 
ume V, BEC occurs in the single particle state $o = iVV, 
having zero momentum. The field operator can be decom- 
posed then in the form $ (r) = y N/V + $' (r). By treating 
the operator $' (r) as a small perturbation, one can develop 
the first or der theory for the excitations of the int eracting 
Bose gases l|Dalfovo et al J 1 19991 ; iBarcelo et alj|200ll ). 

In the general case of a non-uniform and time- 
dependent configuration the field operator in the Heisenberg 
representation is given by $ (f, t) = ip (f, t) + $' (r, t), where 
ip (r,t), also called the condensate wave function, is the ex- 
pectation value of the field operator, tp(r,i) = (<& (r, t)Y It 
is a classical field, and its absolute value fixes the number 
density of the condensate through p(r,t) = \ip (r, t)\ 2 . The 
normalization condition is N = J p(r, t) d 3 f, where N is the 
total number of particles in the condensate. 

The equation of motion for the condensate wave func- 
tion is given by the Heisenberg equation corresponding to 
the many-body Hamiltonian given by Eq. {T}, 

ihl&^t) = [#,#] = 



at 

2m 



V 2 + Vrot (r) + Vext (r) + 



J dr§ + (r, t)V (r -r)$(r, t) 



$ (r, t) . 



(2) 



The zeroth-order approximation to the Heisenberg 
equation is obtained by replacing $ (f, t) with the conden- 
sate wave function ip. In the integral containing the particle- 
particle interaction V (f* — r) this replacement is in general 
a poor approximation for short distances. However, in a di- 
lute and cold gas, only binary collisions at low energy are 
relevant, and these collisions are characterized by a single 
parameter, the s-wave scattering length l a , independently 
of the details of the two-body potential. Therefore, one can 
replace V (r* — r) with an effective interaction V (r* — r) = 
X5 (f 1 — r), where the coupling constant A is related to the 



scattering length l a through A = 47rft 2 Z a /m, where m is the 
mass of the condensed particles. With the use of the effec- 
tive potential the integral in the bracket of Eq. gives 
A \ip (r, t)\ 2 , and the resulting equation is the Schrodinge r 
equation with a quartic nonlinear term (|Chen et al.| [2005). 
However, in order to obtain a more general description of 
the Bose-Einstein condensate stars, we shall assume an ar- 
bitrary non-linear term g (jtp (r, t)| 2 ) = g (p) (jBarcelo et al.l 
l200lf ). 

Therefore the generalized Gross-Pitaevskii equation de- 
scribing a gravitationally trapped rotating Bose-Einstein 
condensate is given by 



2m 



V 2 + Vrot (r) + V Bxt (r) + 



g'M(r,t)\' 



ip{r,t), 



(3) 



where we denoted g' = dg/dp. As for V ex t{r), we assume 
that it is the gravitational potential V, V ex t = V, and it 
satisfies the Poisson equation 



vV 



4-irGpr, 



(4) 



where p m = mp = m\ijj (f,t)\ is the mass density inside the 
Bose-Einstein condensate. 



2.2 The hydrodynamical representation 

The physical properties of a Bose-Einstein condensate de- 
scribed by the generalized Gross-Pitaevskii equation given 
by Eq. (|3]) can be understood much easily by using the 
so-called Madelung representation of the wave function 
ijDalfovo et al.| [l999) . which consist in writing ip in the form 



V»(f,t) = y/p{?,t) 



exp 



(5) 



where the function S (r, i) has the dimensions of an action. 
By substituting the above expression of ip (r, t) into Eq. ©, 
it decouples into a system of two differential equations for 
the real functions p m and v, given by 

dp 



dt 



+ V ■ (p m v) = 0, 



Pm 



di 



-VP 



p m V ( ^) - p m \7 ( ysSL] - X7V Q , 
where we have introduced the quantum potential 

n 2 v 2 A 



/Pm 



V Q = - h — 
2m ^fpm 

and the velocity of the quantum fluid 

m ' 

respectively, and we have denoted 

p ( £hl) = g ' ( Em.) Bn±-g (£™ \ . 

V m J V m / m \m J 

From its definition it follows that the velocity field is ir- 
rotational, satisfying the condition Vx« = 0. Therefore the 
equations of motion of the gravitational ideal Bose-Einstein 



(6) 



(7) 



(8) 



(9) 



(10) 
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condensate take the form of the equation of continuity and of 
the hydrodynamic Euler equations. The Bose-Einstein grav- 
itational condensate can be described as a gas whose density 
and pressure are relate d by a barotropic equation of state 
l|Pethick fc Smitrj|2008h . The explicit form of this equation 
depends on the form of the non-linearity term g. 

When the number of particles in the gravitationally 
bounded Bose-Einstein condensate becomes large enough, 
the quantum pressure term makes a significant contribution 
only near the boundary of the condensate. Hence it is much 
smaller than the non-linear interaction term. Thus the quan- 
tum stress term in the equation of motion of the condensate 
can be neglected. This is the Thomas-Fermi approximation, 
which has been exten sively used for the study of the Bose - 
Einstein condensates (jDalfovo et al.ll 19991 : Ichen et al.ll2005h . 
As the number of particles in the condensate becomes infi- 
nite, the Thomas-Fermi approximation becomes exact. This 
approximation also corresponds to the classical limit of the 
theory (it corresponds to neglecting all terms with powers of 
h, or, equivalently, to the regime of strong repulsive interac- 
tions among particles). From a mathematical point of view, 
the Thomas-Fermi approximation corresponds to neglecting 
in the equation of motion all terms containing Vp and VS. 

In the standard approach to the Bose-Einstein conden- 
sates, the non-linearity term g is given by 



9(Pm) = — \1p\ = —Pm, 



(11) 



wher e u = 4irh 2 l a /m (|Dalfovo et al.l 1 19991 : Ich en et al.1 
2005). The corresponding equation of state of the conden- 
sate is 



P(pm) = Uop: 



m ? 



(12) 



with 



U 



2nh 2 l a 



1.232 x 10° 



Therefore the equation of state of the Bose-Einstein 
condensate with quartic non-linearity is a polytrope with 
inde x n = 1. However, in the case of low dimensional sys- 
tems iKolomeiskv et al l (2000) have shown that in many ex- 
perimentally interesting cases the nonlinearity will be cubic, 
or even logarithmic, in p m . The strong interaction assump- 
tion is valid only if the interaction energy per particle is 
much bigger than the ground-state energy (due to the zero- 
point motion) per particle. This is the case for condensates 
in the dilute limit below two dimensions. But as space di- 
mensionality decreases, it becomes increasingly harder for 
the repulsive particles to avoid collisions. Thus the correla- 
tions between particle dominate, and the quartic nonlinear- 
ity should be replaced b y a more general, power-law term 
|Kolomeiskv et alj[200oh . Hence more general models, with 
the non-linearity term of the form g (p m ) = aPm> where a = 
constant and T = constant, can also be considered. In this 
case the equation of state of the gravitational Bose-Einstein 
condensate is the standard polytropic equation of state, 
P (pm) = a(r — 1) p^, and the structure of the static gravi- 
tationally bounded Bose-Einstein condensate is described by 
the Lane-Emden equation, (l/£ 2 ) d (£ 2 d6>/d£) /d£ + 6 n = 0, 
where n = 1/ (r — 1) and 8 is a dimensionless variable de- 
fined by p = Pcm0 n , where p cm is the central density of 



the condensate. The dimensionless radial coordinate £ is 

I 1/2 



(n + l)Kpl£ 



V47TG 



defined by the relation r — 

Hence Bose-Einstein condensate dark matter can generally 
be described as fluid satisfying a polytropic equation of state 
of index n. 

In the following we will consider only the case of the 
condensate with quartic non-linearity. In this case the phys- 
ical properties of the condensate are relatively well known 
from laboratory experiments, and its properties can be de- 
scribed in terms of only two free parameters, the mass m 
of the condensate particle, and the scattering length l a , re- 
spectively. 



2.3 Dark matter as a Bose-Einstein condensate 

In the case of a static Bose-Einstein condensate, all physical 
quantities are independent of time. Moreover, in the first ap- 
proximation we can also neglect the rotation of the conden- 
sate, by taking V ro t = 0. Therefore the equations describing 
the static Bose-Einstein condensate in a gravitational field 
with potential V take the form 



VV = -i-nGpsE- 



(14) 
(15) 



These equations must be integrated together with the 
equation of state P = P (pbe) = Uop% E , and some appro- 
priately chosen boundary conditions. The density distribu- 
tion pBE of the static gravitationally bounded single com- 
ponent dark matter Bos e-Einstein condensate is given by 
ijBoehmer fe Harkd 120071 ) 



Pbe (r) 



( c ) sin kr 

pBE ~kV 



(16) 



where k — yj Gm 3 /h 2 l a and p EE is the central density of 



the condensate, p EE = pse(0). The mass profile niBE(r) = 
4-7T J Q r pBE(r)r 2 dr of the Bose-Einstein condensate galactic 
halo is 



4 (c) 
m BE (r) = ^gPr 



' sin kr 
. kr 



— cos kr 



(17) 



with a boundary radius Rbe- At the boundary of the dark 
matter distribution Pbe(Rbe) ~ 0, giving the condition 
kRsE = 7r, which fixes the radius of the condensate dark 
matter halo as Rbe = ~K\jTi 2 l a /Grrv i . The tangential veloc- 
ity of a test particl e moving in the condens ed dark halo can 
be represented as (jBoehmer fc Harkoll2007h 



vIe (r) = 



Gm,BE(r) 4arGpg E / sin kr 



k 2 



kr 



cos kr 



(18) 



The mass of the particle in the condensate can be ob- 
tained from the radius o f the dark matter halo in the form 
ijBoehmer fc Harkdl2007h 

V GR BE J 

2.58 x 10~ 30 [l a (cm)] 1/3 [R BE (kpc)]- 2/3 g » 
6.73 x 10" 2 [l a (fm)] 1/3 [Rbe (kpc)]" 2/3 eV. (19) 



Cosmological perturbations in Bose-Einstein condensates 5 



From this equation it follows that m is of the order of 
eV. For l a « 1 fm and Rbe ~ 10 kpc, the mass is of the or- 
der of m ~ 14 meV. For values of l a of the order of l a ~ 10 6 
fm, corresponding to the values of l a observed in terrestrial 
laboratory experiments, m w 1.44 eV. These values are per- 
fectly consistent with the limit m < 1.87 eV obtained for 
the mass of the condensate partic le from cosmological con- 
siderations (Fuk uvama et al1l2008l ). 



3 POST-NEWTONIAN HYDRODYNAMICS 
OF THE BOSE-EINSTEIN CONDENSATES 

In order to study gravitational effects on the evolution of 
Bose-Einstein condensate dark halos a full general relativis- 
tic treatment is needed. The equations of motion of the con- 
densate are obtained from the conservation of the energy- 
momentum tensor, T£„ = 0, with ; denoting the covariant 



derivative with respect to the metric g^, and 
TS = ( Pm c 2 + P) u v u^-P8Z, 



(20) 



where is the four-velocity of the fluid, satisfying the con- 
ditions u^u^ — 1 and M^it M = 0, respectively. By taking the 
covariant divergence of Tjf we obtain the equation 



(pmC 2 + P) UvU^ + [p m C + P) M„ |( 
(pmC 2 + P) U V U^ = P >u , 



(21) 



where a comma denotes the ordinary derivative with respect 
to the coordinate . By contracting Eq. (pTjl with u", we 
obtain 



J +P) u»=P„u. 



(22) 



In the Newtonian limit of small condensate velocities 
the four- velocity is given by u M = (l,v/c), where v is the 
three-velocity of the condensate. The four-divergence of the 
four- velocity is given by = {l/y/—g) d (yf—gu 11 ) /dx^, 
where (— g) is the determinant of the metric tensor. In the 
Newtonian limit we assume that ^J—g — > 1, that is, the devi- 
ations from the Minkowski type geometry are small. Under 
these assumptions, from Eq. (|22[) we obtain the equation of 
continuity of the Bose-Einstein condensate as 



+ v,. <*.*)+ £v r 



v = 0, 



(23) 



where all differential operations are considered with respect 
to the physical coordinate f. By contracting Eq. (|2ip with 
the projection operator h v a — 8^ — u a u v , with the property 
0, we obtain the relativistic Euler equation of mo- 



1 iuii as 

2 



(p m c 2 + P) U q;m U M = P t0l — P, v UU a . 



(24) 



In the Newtonian approximation the generalized Euler 
equation of motion becomes 



-VrV 



cVrP + Pv 
PmC 2 + P 



(25) 



The gravitational potential V satisfies the generalized 
Poisson equation, 



4tvG (p m + 3 



(26) 



Eqs. (|23|) . (|25|) and (|26|) represent the basic equations 
describing the dynamics of a gravitationally bounded Bose- 
Einst ein condensate in the first Post - Newtonian approxima- 
tion (|McCrealll95ll; lHarrisonlll965l; iLima et alj 1 19971 ; iReisI 



l2003l ; 1 Abramo et al.H2007l ; iPace et al.ll2010l ) 



4 COSMOLOGICAL DYNAMICS OF 
BOSE-EINSTEIN CONDENSATES 

The Bose-Einstein condensation takes place when particles 
(bosons) become correlated with each other. This happens 
when their wavelengths overlap, that is, the thermal wave- 
length At = y 2ivh 2 /raksT is greater than the mean inter- 
particles distance a, At > a. The critical temperature for 
the condensation to take place is T cr < 2nh 2 n 2 ' 3 /mkB 
l|Dalfovo et alj [l999). On the other hand, cosmic evolution 
has the same temperature dependence, since in an adiabatic 
expansion process the densi ty of a matter dominated Un i- 
verse evolves as p oc y 3 / 2 ||Fukuvama fc Morikawall2009l ). 
Therefore, if the boson temperature is equal, for example, 
to the radiation temperature at z = 1000, the critical tem- 
perature for the Bose-Einstein condensation is at present 
T cr = 0.0027^ l|Fukuvama fc Morikawal |2009| ). Since the 
matter temperature T m varies as T m oc a~ 2 , where a is 
the scale factor of the Universe, it follows that during an 
adiabatic evolution the ratio of the photon temperature P 7 
and of the matter temperature evolves as T 7 /T m oc a. Us- 
ing for the present day energy density of the Universe the 
value p^ = 9.44 x 10~ 30 g/cm 3 , BEC takes place provided 
that the boson mass s atisfies the restriction m < 1.87 eV 
IjFukuvama et al.ll2008l ). Thus, once the temperature T cr of 
the boson is less than the critical temperature, BEC can 
always take place at some moment during the cosmologi- 
cal evolution of the Universe. On the other hand, we expect 
that the Universe is al ways under critical temperatur e, if it 
is at the present time IjFukuvama fc Morikawal 120091 . An- 
other cosmological bound on the mass of the condensate 
particle can be obtained as m < 2.696 (g d /g) {T d /T cr f eV 
IjBovanovski et alj |2008). where g is the number of internal 
degrees of freedom of the particle before decoupling, g d is 
the number of internal degrees of freedom of the particle at 
the decoupling, and T d is the decoupling temperature. In the 
Bose condensed case T d /T c < 1, and it follows that the BEC 
particle should be light, unless it decouples very early on, at 
high temperature and with a large g d . Therefore, depend- 
ing on the relation between the critical and the decoupling 
temperatures, in order for a BEC light relic to act as cold 
dark matter, the decoupling scale must b e higher than the 
electroweak scale l|Bovanovski et al]|2008l ). 

The set of equations Eqs. (|23p . (|25|l and (|26p admits a 
homogeneous and isotropic cosmological background solu- 
tion with p m — pb(t) and P = p, (t). In this case the fluid's 
velocity is given by 

vt = V, (27) 

and the evolution of the scale factor a is determined by the 
Friedmann equations, 

(28) 



3^- = 3H 2 = 8ivGp b 
and 
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4ttG 
~3~ 



Pb - 



Pb 



(29) 



respectively, where we have denoted H — a/a. The continu- 
ity equation Eq. ()23[) reduces to 



= 0. 



(30) 



In the case of the Bose-Einstein condensates the equa- 



tion of state is given by Pb = Uop b , and Eq. 
integrated immediately to obtain 



Pb (a) 



C 



CU /c 2 



3011 can be 



(31) 



where C is an arbitrary constant of integration. By assum- 
ing that the present-day density of the Bose-Einstein con- 
densate, p m ,o, is obtained for a value a — ao of the scale 
factor, we obtain C = p m ,o«o/ (l + pm,oUo/c 2 ^, and the 
background cosmological density of the condensate can be 
written as 



pb (a) = 



Po 



Uo (a/a )' - po ' 
where we have denoted 

p m ,oUo/c 2 



Po = 



1 + pmfiUo/c 2 



(32) 



(33) 



The energy density of the Bose-Einstein condensate di- 

1/3 

p . The equation determining the time evo- 



verges as a 
lution of the scale factor is given by 



— = Ho\/Qbe 
at 



\J (a/a ) 3 - po 
where we have denoted 

Qbe,o "f 



Q 



l + p m ,oU /c 2 1 + O.BE,OPcr,oU /c 2 



(34) 



(35) 



where Hbe.o = Pm,o/Pcr,o is the present day density param- 
eter of the Bose-Einstein condensate, p cr ,o = 3Hq/8ttG is 
the present day critical density of the Universe, and Ho is 
the present day value of the Hubble parameter, respectively. 

Eq. (|34p can be integrated immediately to give the time 
evolution of the scale factor of the Bose-Einstein condensate 



VHeHHo (t-Ci) 



(-) 

\aoJ 



Po 



- yfpo arctan 



(a /ao) 



Po 



po 



(36) 



where Ci is an arbitrary constant of integration. The con- 
stant Ci can be determined from the condition t = when 
(a/ao) 3 = po, thus obtaining C\ = 0. Therefore the time 
evolution of the scale factor is described by the equation 



t 



3\ / TTbe 
yfpo arctan 



\J {a/ao) 



Po 



(a/ao) 



Po 



Po 



(37) 



where we denoted tu = 1 /Ho. 

In the case of the standard dark matter models, dark 




t/'H 

Figure 1. Time evolution of the scale factor of the pressureless 
dark matter, described by the Einstein-de Sitter metric (solid 
curve), and of the Bose-Einstein condensate dark matter, for dif- 
ferent values of po: po = 10~ 6 (dotted curve), po = 10~ 3 (dashed 
curve), po = 0.005 (long-dashed curve), and po = 0.01 (ultra-long 
dashed curve). 



matter is assumed to be a pressureless fluid, and the back- 
ground cosmological evolution is described by the Einstein- 
de Sitter model, with the scale factor given by a/ao = 
(9f2rjA/,o/4) 1//3 (t/tn) 2 ^ 3 , where Qdm,o is the present day 
density parameter of the dark matter, and we have assumed 
that a(0) = 0. In the following for the Hubble constant we 
adopt the value Hp = 70 km/s/Mpc = 2.273 x 10" 18 s" 1 
l|Hinshaw et alj|2009h . giving for the critical density a value 
of pc,o = 9.248 x 10"" 
resented as 



/cm . The constant po can be rep- 



po 



1.266 x fl B E,o x (m/1 meV)" 3 x (/ a /10 9 fm) 



1 + 1.266 x Q BE ,o x (m/1 meV) -3 x (Z a /10 9 fm) 
while for Qbe we obtain 



(38) 



Of 



Q 



BE,0 



1 + 1.266 x Q.be,o x (m/1 meVy ' x (Z a /10 9 fm) 



.(39) 



The two parameters po and Q.be, describing the global 
properties of the condensate, are related by the relation 
Q.be — Odm,o x (1 — po). 

By assuming that the entire existing dark matter is 
in the form of a Bose-Einst ein condensate, it fol lows that 
Mbe,o w n DM ,a w 0.228 (|Hinshaw et alj 120091 ). By as- 
suming that m — 1 meV and l a = 10 10 fm, we obtain 
po = 0.7426, while £1 B e = 0.05866. For these values of the 
physical parameters of the condensate the energy density 
of the dark matter diverges for a/ao — > 0.9056. The time 
evolution of the scale factor a for the Bose-Einstein con- 
densate dark matter, for different values of the parameters 
m and l a , and for the pressureless dark matter, are repre- 
sented in Fig. [TJ respectively. The cosmological dynamics of 
the condensate dark matter shows significant differences as 
compared to the standard pressureless dark matter model, 
with the condensate expanding much faster than the cosmo- 
logical fluid of the standard ACDM model, with the speed 
of expansion increasing with increasing po. 

In the case of a Universe filled with dark energy, radi- 
ation, baryonic matter with negligible pressure, and Bose- 
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Figure 2. Time evolution of the scale factor of a Universe 
filled with dark energy, baryonic matter and with Bose-Einstein 
condensate dark matter, respectively, for different values of po: 



po = 10 (dotted curve), po 



10" 



(dashed curve), po = 0.05 



(long-dashed curve), and po = 0.1 (ultra-long dashed curve). The 
solid curve represents the time evolution of the scale factor in the 
standard ACDM model, in which the dark matter is pressureless. 



Einstein condensed dark matter, respectively, the time evo- 
lution of the scale factor is given by 



fir 



Ida = H 

a dt °V ( a /a ) 3 ' (a/a ) 4 ' (a/a f 



Po 



+ n A , (40) 



where Qb,o, firad.o, and Q\ are the present day values of the 
density parameters of the baryonic matter, radiation, and 
dark energy, respectively. For SIb,o, ^rad,o, an d £l A we adopt 
the numerical val ues SIb.o = 0.0456, i l r ad,o = 8.24 x 10~ 5 
and tt A = 0.726 ijHinshaw et all 120091 ). In the case of the 
standard ACDM cosmological, model the Friedmann equa- 
tion describing the evolution of the Universe containing 
baryons, pressureless dark matter, radiation, and dark en- 
ergy, is given by 



1 da 
a dt 



H 



(a/ao) 3 



(a/a ) 



(41) 



The time evolutions of the scale factors for Universes 
containing BEC dark matter and standard pressureless dark 
matter are represented, for different values of the BEC pa- 
rameter po, in Fig. [2] The presence of the condensate dark 
matter changes the global cosmological dynamics of the Uni- 
verse, and the magnitude of the changes increases with the 
increase of the BEC parameter pa- 



in Eqs. (|42[l - (|45fl the index b denotes the background 
quantities. Substituting these equations into the continuity 
equation Eq. (123[l we obtain 



+ ( P6 + 50 Vr-U + Vr-(S P V b ) + ^Vr-ti b =0.(46) 

The variation of the equation of motion Eq. (|25p gives 
dit 



dt 



—V r <p — 



+ (Vb ■ Vr) U + (u ■ Vr) Vb 
V r 5P + P b u/C 2 



Pb + Pb/c 2 



(47) 



The Poisson equations for the perturbation of the grav- 
itational potential, obtained by perturbing Eq. I]26[l. can be 
written as 

(48) 



In order to describe the cosmological evolution we make 
a change to the comoving coordinate system, so that r = aq, 



V 9 = V = aV r , and 
f d\ d ( d \ a „ 



(49) 



respectively. To simplify the notation we define the param- 
eters w = Pb/pbC 2 and c?** = 5P/8p, respectively, which 
generally are functions of the time only. We also introduce 
the density contrast as 8 = Sp/pb- The time derivative of 
the background pressure is related to the speed of sound of 
the background condensate c 2 = dPb/dpb by the relation 



Pb = 



-3 r^j c 2 p b (i + w). 



(50) 



Therefore the perturbation equations Eqs. (|43[) - (|45[1 can 
be written as 

S + 3H (^-w S jS+^-^V-u = 0, (51) 

du /„ „c? \ 1„ c 2 ff 1 1 

— + (1-3-f )-u+-Vv+^f-— — V-i = 0, 52 

at \ cr j a a r nl + » 



V 2 ¥> = AnGa 2 p b ( 1 + 3-^ ] 6. 



In the following we denote 



C eff 



(53) 



(54) 



and 



5 COSMOLOGICAL PERTURBATIONS OF AN 
EXPANDING BOSE-EINSTEIN 
CONDENSATE 

In the gravitationally bounded Bose-Einstein condensate we 
assume small perturbations of the physical quantities around 
the homogeneous background of the form 



Pm (r, t) = p b (t) + 5p (r, t) , 
P(r y t) = P b (t) + 5P(r,t), 
V(r,t) = V b + <p(f,t), 
V (f, t) = v b +u (f, t) . 



(42) 
(43) 
(44) 
(45) 



Q-s = -s — W, 

cr 



(55) 



respectively. a s is related to the time derivative of w by the 
relation 



= -3Ha„ 



(56) 



l + w 

By taking the time derivative of Eq. (|51|l , the divergence 
of Eq. (J52J, by eliminating V-w by using the perturbed equa- 
tion of continuity, and with the use of Eq. (|53[) . we obtain 
the equation giving the evolution of the density contrast as 



A -Ml j -' +l)*+l"' 2 x 
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9w — 2w 



. cl ff 2 d 
2(l + 6w )^ + --a eff 



- 1 



S = 



-eff 



-AS. 



(57) 



and for its derivative for BEC dark matter requires the full 
investigation of the cosmological dynamics from the reheat- 
ing era, by taking into account the dark matter condensation 
and pressure effects. 



Changing the independent variable from the time t to 
the scale factor a using the relations d/dt = aH(a)d/da and 
d 2 /dt 2 = a 2 H 2 d 2 /da 2 - [(1 + 3w)aH 2 /2] 8/ da, we obtain 
the evolution equation in the form 



2 d*5 
a — — + 3a 

oa 2 



9w — 2w 



5 1 
— w H — 
2 2 



85 
da 



c 2 tt 

2 (1 + 6w) -^f + 2a- 

cr < 



a 2 H 2 



-AS. 



(58) 



Eq. (|57[) is different from the perturbation equation 
obtained in the N ewto n ian c osmol ogy with pressure b y 
iLima et all (|l997l ). iReisI (|2003l ), and lAbramo et alj (|2007l ). 
respectively. The reason is that we have included in our anal- 
ysis the term Pbu/c 2 , which was neglected in the previous 
studies. This term generates the new term — (3c 2 /c 2 ) Hu 
in the left hand side of Eq. (|52|) . which modifies the final 
perturbation equation. On the other hand, in the present 
approach the term 5P (f,t) /c 2 was neglected. 

In order to numerically integrate Eq. (|57[) or Eq. (|58|) 
we have to chose some physically appropriate initial condi- 
tions. In the current standard model for structure forma- 
tion in the Universe, it is supposed that quantum fluctu- 
ations were generated during an initial period of inflation. 
These fluctuations inflated up to super-horizon scales, pro- 
ducing a near scale-invariant, and near Gaussian, set of pri- 
mordial potential fluctuations. At the end of inflation, the 
Universe is reheated, and particles and radiation are pro- 
duced. In this hot ear ly phase, cold thermal relics (dark mat- 
ter) a re also formed |Peebles fc Ratrall2003l : TPadmanabhanl 
2003). Dark matter particles interact gravitationally, and 
possibly through the weak interaction. Therefore in order 
to obtain some physically realistic initial conditions for the 
state of the Universe during large scale structure forma- 
tion, one needs to evolve cosmological perturbations, start- 
ing from initial conditions, deep inside the radiation epoch, 
and far outside the Hubble radius. Initial conditions for pho- 
tons, neutrinos, cold dark matter and baryons have been ob- 
tained, in the framework of the standard ACDM cosmolog- 
ica l model, in both the synch ronous and Newtonian gauges, 
bv lMafc Bertschingeri l|l995h . In the conventional method, 
the power spectrum of the matter fluctuations in the Uni- 
verse is computed by numerically solving the Boltzmann 
equation. The power spectrum is usually obtained in the 
linear theory, and then extrapolated to the present epoch. 
In the standard cosmology, in first-order Eulerian perturba- 
tion theory, all modes evolve independently, and the power 
spectrum can be scaled back to the initia l epoch via the 
growth function |Ma fc Bertschingeri fl9 95f). Moreover, the 
effect of the dark matter pressure is generally ignored in the 
conventional methods of generating initial conditions. On 
the other hand, the exact moment in the history of the Uni- 
verse when the Bose-Einstein condensation occurred is not 
known. That's why obtaining the rigorous and physically 
well motivated initial conditions for the density contrast 8 



6 COSMOLOGICAL EVOLUTION OF SMALL 
PERTURBATIONS IN A BOSE-EINSTEIN 
CONDENSATE 

By taking into account the equation of state of the Bose- 
Einstein condensate, we immediately obtain w — Pb/pbC 2 — 
/c 2 ) p b = po/ \(a/a ) 3 - pol , and c 2 s /c 2 = c\ ss /c 2 



(Uo/c 2 ) p b = p / [(a/a ) 3 - po] , and c 2 s /c 2 = c 2 eff /c 2 = 2w, 
respectively. The conditions c 2 /c 2 ^ 1 and c^jj/c 2 ^ 1 im- 
poses the constraint (a/ao) 3 ^ 3po, and in the following we 
will consider that the model considered in the present pa- 
per is valid only for this range of values of the scale factor. 
Hence for the time evolution of the linear density perturba- 
tions of the Bose-Einstein condensate dark matter we obtain 
successively 

2 d 2 5 3 ,„ 



§( 



1 + 6w + 15w 



dS_ 

da 

2a^)5- 



da J 



0. 



(59) 



and 



, d 2 5 
da 2 



1 + 



2 

6po 



d5_ 

da 



+ 



15p 2 



+ 



6poa 3 



PO 



(a 3 - po) (c 



po) 



5 = 0, (60) 



respectively, where we have denoted a = a/ao. In the limit of 
large a, when (a/ao) 3 >> po, w 

2 d 2 5 3 dS 3 X n 
Q rf^ + 2 a ^-2^ ' 



with the solution 

J(a) „ Cl (-2-)+ <*(-•)" 

V ao / V ao / 



3/2 



0, and Eq. (|60[) becomes 
(61) 

(62) 



where C\ and C2 are arbitrary constants of integration. 
Hence in the limit of a pressureless fluid we recover the stan- 
dard general relativistic result. In the limit of small a, so 
that a 3 — > 3po, we can approximate Eq. (|60[) as 



2 d 2 5 3 dd 
a — — ■ + -a- 

dar 4 da 

with the solution 



147 , 



0, 



5 (a) w Ci 



(-) 

\a J 



(l+vTL77)/8 



+ c 2 



(63) 



(64) 



where C 1 and C 2 are two arbitrary constants of integration. 
By introducing w = po/ (a 3 — po) as a new indepen- 



■ \W A (1 + w) 2 ^Ar + ~w (1 + w) (1 + 3w) ^~ 
dw 2 2 dw 



dent variable, Eq. (|60|l can be written as 



2 , ^2 d 2 5 

t 

i (l + 12w + 21w 2 ) 



6 = 0. 



(65) 



By representing the density contrast S as 8(w) = w~ (1 + 
u(w), it follows that the new function u(w) satisfies 
the equation 



-5/3 
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d 2 u 



du 



Gw(l+w) — — + (1 - 9w) 15u = 0. 

aw 2 dw 



(66) 



Therefore the general solution of Eq. (|65[) can be obtained 

as 



8(w) = w~ 1/3 (l + w) 



-5/3 



C[ 2 Fi 



5 + V6E + 1 
r 6 ; 



r,i 5/6 rp 
C 2 w 2 Fx 



_5_ 

12 



4 

65 5 



12 + 4 



65 11 
6 



(67) 



where 2F1 (a, b; c; z) is the hypergeometric function, 
2 F 1 (a,b-c;z) = J^Zo (b)^V(c)fcfc!, | z | < 1, and CJ 
and C 2 are arbitrary constants of integration. The constants 
of integration can be determined from the initial conditions. 
When a 3 = 3p , w = 1/2, and w(l/2) = (3 5/3 /4) Si and 
u'(l/2) = (3 5/3 /4) [5'i + 165 4 /9], where we have denoted 
<5; = 5(1/2) and 8[ — S'(l/2), respectively. In the limit of 
large a, w — > 0, and the density contrast can be approxi- 
mated as 

8(w) w 0.00231255 (347.503^ - 250.317^) w~ 1/3 + 

5C 



0.00925 



1533.98 



j I - 2366.515i 



0.0308341 (250.317(5- - 347.503(5;) w 2/3 + 0.00925 x 



5019.; 



1533.98 



1.66667 x 



y ) - 2366.515; 



+ 



3253.9 0.888(5,: + 



5 A 

{5094.725; - 1.66667 x 
2971.925; - 1877.38 | 

2.22222 

3218.36 



w 3/2 + 0.00925022 x 



5'i 

1 + 2 



198.1285; - 125.159 



+ 



0.8885, + j w5/3 + ( w7/3 ) . w -»• 0- 



(68) 

Near the initial state w = 1/2, the density contrast can be 
approximated as 



(5 (to) 



+ 2.39916 x 10~ u x 
2 



(7564365; - 5789055;) - 
2.97184 (-1.102725, + 0.864771$) ( to - -J + 
,u;-H/2. (69) 



Since Eq. (|60|l is valid only in the linear regime of small 
perturbations, we assume that the initial value of the per- 
turbation, 8 (ai/ao), occurring for a value a = 04 of the 
scale factor, satisfies the condition 5 (ai/ao) « 1. Since the 
equation describing the perturbations is a second order dif- 
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Figure 3. Evolution of the density perturbations 8 (a/ao) as a 
function of a/ao of the pressureless dark matter, described by the 
Einstein-de Sitter metric (solid curve), and of the Bose-Einstein 
condensate dark matter, for different values of po: po = 10~ 3 
(dotted curve), po = 2 X 10~ 3 (dashed curve), po = 3 X 10~ 3 
(long-dashed curve), and po = 5 X 10 -3 (ultra-long dashed curve). 
In all cases the initial conditions are «5(3p ) 1/3 = 10" 5 and 
8' (3po) 1/3 = 10~ 5 , respectively. 
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Figure 4. Evolution of the density perturbations <5(a/ao) as a 
function of a/ao °f the pressureless dark matter, described by the 
Einstein-de Sitter metric (solid curve), and of the Bose-Einstein 
condensate dark matter, for different values of po: po = 10 -3 
(dotted curve), po = 2 X 10~ 3 (dashed curve), po = 3 X 10~ 3 
(long-dashed curve), and po = 5 X 10 — 3 (ultra- long dashed curve). 
In all cases the initial conditions are «5(3p ) 1/3 = lCT 5 and 
8' (3po) 1//3 = 1.5, respectively. 



ferential equation, two initial values have to be given, one 
for the initial perturbation 5 (ai/ao), and one for the initial 
rate of evolution of the perturbation, 8' (ai/ao). We consider 
two cases, namely, the case of a perturbation with an ini- 
tial low evolution rate, of the order of &' (ai/ao) = 10~ 5 , 
and the case of a perturbation with a very high initial evo- 
lution rate, 8' (ai/ao) = 1.5, respectively. The comparison 
between the evolution of the evolution of the linear pertur- 
bations for pressureless dark matter in an expanding Ein- 
stein - de Sitter cosmological background, and the evolution 
of the density perturbations in a Bose-Einstein condensate 
dark matter dominated Universe is presented in Figs. [3] and 
U respectively. 

As one can see from Figs. [3] and [4] for all initial con- 
ditions, in the case of the Bose-Einstein, for a given a, the 
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amplitude of the density contrast is higher as compared to 
the case of the standard dark matter model. The condensate 
enters more rapidly in the non- linear phase (S » 1) than 
the pressureless dark matter. Thus the presence of the Bose- 
Einstein condensate dark matter can significantly accelerate 
the process of cosmic structure formation. 



7 DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the global cosmo- 
logical evolution and the evolution of the small cosmologi- 
cal perturbation in a Bose- Einstein dark matter condensate. 
The basic equation describing the evolution of the small per- 
turbations in the Post-Newtonian regime was obtained, and 
its solutions have been studied by using both analytical and 
numerical methods. The evolution of the density perturba- 
tions of the condensate has been compared to the evolution 
of the small cosmological perturbations in a pressureless 
fluid evolving in an Einstein-de Sitter cosmological back- 
ground. Depending on the numerical values of the physical 
parameters describing the condensate (the mass of the par- 
ticle and the scattering length, respectively), significant dif- 
ferences could appear in the evolution of the Bose-Einstein 
condensate dark matter halos, as compared to the standard 
pressureless dark matter models. These differences appear 
at both the level of the global cosmological evolution, and 
of the behavior of the small perturbations in the dark mat- 
ter fluid, and they could have fundamental implications for 
the formation of the large-scale structure in the Universe. 

One of the most important problems present day cos- 
mology faces is the problem of the galaxy formation. To 
explain galaxy formation the evolution of the linear density 
and temperature perturbations in a Universe with dark mat- 
ter, baryons, and radiation must be computed. For pressure- 
less dark matter the evolution of the perturbations of all cos- 
mic components, from cosmic reco mbination until th e epoc h 
of the first galaxies, was obtained in iNaoz fc Barkanal|2005l) . 
The evolution of sub-horizon linear perturbations can be de- 
scribed by two coupled, second-order differential equations, 
with the pressureless dark matter interacting gravitation- 
ally with itself, and with the baryons, while the baryons 
experience both gravity and pressure. Starting from very 
low values on sub-horizon scales, the baryon density pertur- 
bations gradually approach those in the dark matter, and 
the temperature perturbations approach the value expected 
for an adiabatic gas. The presence of the baryons does not 
modify significantly the evolution of the dark matter per- 
turbations. By including the effect of the BEC pressure in 
the perturbations equations for baryons and dark matter, a 
more general (and realistic) description of the galaxy forma- 
tion process can be obtained. The presence of the BEC mod- 
ifies the dynamical evolution of the baryons, and the growth 
of linear perturbations, which provide the initial conditions 
for the formation of galaxies. In the BEC condensate model 
the dark matter perturbations grow more rapidly than in 
the standard cosmology, and therefore this could lead to a 
much faster growth rate of the baryonic perturbations, ac- 
celerating the galaxy formation process. 

A major recent experimental advance in the study of 
the Bose-Einstein condensation processes was the observa- 
tion of the collapse and subsequent explosion of the conden- 



sates l|Rvbin et alj|2003 ) . A dynamical study of an attractive 
85 Rb BEC in an axially symmetric trap was done, where 
the interatomic interaction was manipulated by changing 
the external magnetic field, thus exploiting a nearby Fesh- 
bach resonance. In the vicinity of a Feshbach resonance the 
atomic scattering length a can be varied over a huge range, 
by adjusting an external magnetic field. Consequently, the 
sign of the scattering length is changed, thus transforming a 
repulsive condensate of 85 Rb atoms into an attractive one, 
which naturally evolves into a collapsing and exploding con- 
densate. From a simple physical point of view the collapse of 
the Bose-Einstein condensates can be described as follows. 
When the number of particles becomes sufficiently large, so 
that N > N c , where N c is a critical number, the attrac- 
tive inter-particle energy overcomes the quantum pressure, 
and the condensate implodes. In the course of the implosion 
stage, the density of particles increases in the small vicin- 
ity of the trap center. When it approaches a certain critical 
value, a fraction of the particles gets expelled. In a time pe- 
riod of an order of few milliseconds, the condensate again 
stabilizes. There are two observable components at the fi- 
nal stage of the collapse: remnant and burst particles. The 
remnant particles are those which remain in the condensate. 
The burst particles have an energy much larger than that of 
the condensed particles. There is also a fraction of particles, 
which is not observable. This fraction is usually referred to 
as the missing particles. 

The scattering length l a is d efined as the zero-e nergy 
limit of the scattering amplitude / l|Dalfovo et al ][l99l). De- 
pending on the spin dependence of the underlying particle 
interaction, the scattering length may in general be also spin 
dependent. The spin independent part of the quantity is re- 
ferred to as the coherent scattering length l a . The scattering 
lengths can be obtained for some systems in the laboratory, 
but for dark matter it is unknown. Another essential param- 
eter is the mass m of the condensate particle, which, due 
to the lack of information about the physical nature of the 
dark matter, is a free parameter, which must be constrained 
by observations. Due to the lack of any physical informa- 
tion about the numerical values of these two fundamental 
parameters, in the numerical estimations performed in the 
present paper we have given different numerical values to a 
combination of these two basic quantities. 

Since Bose-Einstein condensates are less stable with re- 
spect to perturbations than usual non-condensate matter, 
we expect that in such a condensate evolution of the pertur- 
bations and the subsequent collapse could take place much 
faster than in the usual non-condensed matter. This would 
strongly affect the formation of the large scale structure in 
the early Universe. In this paper we have provided some 
basic theoretical tools necessary for the in depth compari- 
son of the predictions of the condensate model and of the 
observational results. 
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